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Abstract 

Locally any completely integrable system is maximally superintegrable system such as 
we have the necessary number of the action-angle variables. The main problem is the 
construction of the single- valued additional integrals of motion on the whole phase space 
by using these multi-valued action-angle variables. Some constructions of the additional 
integrals of motion for the Stackel systems and for the integrable systems related with 
two different quadratic r-matrix algebras are discussed. Among these system there are 
the open Heisenberg magnet and the open Toda lattices associated with the different root 
systems. 

1 Introduction 

In classical mechanics a Hamiltonian system on a 2n-dimensional phase space M is called 
completely integrable in Liouville's sense if it possesses n functionally independent integrals of 
motion Hi , . . . , Hn in the involution: 

^ = {i/,ffJ-0, {H,,H,}^Q, 

where H = Hi is the Hamilton function and {., .} is the Poisson bracket on M. 

A superintegrable system is a system that is integrable in the Liouville sense and that, in 
addition to this, possesses more functionally independent integrals of motion than degrees of 
freedom. If the number of independent integrals takes the value 2n — 1, then the system is 
called maximally superintegrable [M] . There are three classic and well-known examples of such 
systems, namely, the free particle that can be considered as trivial, the Kepler problem, and 
the harmonic oscillator with rational frequencies. In these cases all the orbits become closed 
for the case of bounded motions. 

The Liouville classical theorem on completely integrable Hamiltonian systems implies that 
almost all points of the manifold M are covered by a system of open toroidal domains with the 
action-angle coordinates / = {Ii, . . . , Ik) and — (wi, . . . , w„) 

{Lj,h}^{LO,,LUk}^0, {/,,L^fe} = %. (1.1) 

The independent integrals of motion Hi , . . . , Hn are functions upon the independent action 
variables /i , . . . , /„ and the corresponding Jacobian does not equal to zero 



detJ=^0, where J„- = ' " ' " ■ (1-2) 



Let us introduce n functions 



'.-E(J")..-'=' (1-3) 



such that 

n 

= ^J,fc (J-i),^. (1.4) 



fe=i 
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The (n — 1) functions 02, ■ • ■ , are integrals of motion 

^ = 0, j^2,...,n, 

which are functionally independent on n functions Hi{I), . . . , Hn{I). So, in classical mechanics 
any completely integrable system is superintegrable system in a neighborhood of any regular 
point of M. It means that the Hamiltonian H — Hi has 2(n — 1) integrals of motion H2, ■ ■ ■ , 
and (f>2, ■ ■ ■ , (t>n on any open toroidal domain. 

Sometimes the action-angle variables are global variables on the whole phase space M and, 
therefore, we have superintegrable systems on M. For instance, the global action-angle variables 
for the open and periodic Toda lattices are discussed in [Sj. 

In generic case the action variables ujk are multi-valued functions on the whole phase space 
M. In this case if we have some algebraic integral of motion K, it has to be some function 
on the action-angle variables. In this paper we discuss some possibilities to get polynomial 
integrals of motion from the multi-valued angle variables. 



2 The Stackel systems. 

The system associated with the name of Stackel [13l HSl [16] is a holonomic system on the phase 
space R^", with the canonical variables q = {qi, . . . ,qn) and p = (pi, . . . ,p„): 

n 

r2 = ^ dp-j A dq-j , {P],qk} = Sjk ■ (2.1) 

The nondegenerate n x n Stackel matrix S, whose j column depend on coordinate qj only, 
defines n functionally independent integrals of motion 

n 

Hk^Y.^S-'^ik{p',+U,{q,)). (2.2) 
From this definition one immediately gets the separated relations 

n 

p] = J2HkSk,{q,)~U,{q,). (2.3) 

fc=l 

It allows reducing the solution of the equations of motion to a problem in algebraic geometry 
[U m [TS] . We can regard each separated equation (|2.3p as equation defining the hyperelliptic 
curve 

Cr- ^i^=P,{X,), (2.4) 

where fij and Xj are functions on pj and qj only and Pj{Xj) are polynomials. 

The functionally independent action variables Ik = Hk (12. 2p have the canonical Poisson 
brackets (jl.ip with the following angle variables 



which consist of integrals of first kind Abelian differentials on the hyperelliptic curves Cj (|2.4p 
[U m [T21 US]. For the coinciding curves we have to use different A-cycles on the Riemann 
surfaces in the sums (|2.5p . 

In generic case the action variables (|2.5p are the sum of the multi- valued functions "dij. 
However, if we are able to apply some addition theorem to the calculation of uji (|2.5p 



t^. = ^*y(Pj,9j) = e,(if,), (2.6) 
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where Qi is a multi- valued function on the algebraic argument Ki {p, q) then one gets algebraic 
integrals of motion Ki{p, q) such that 

{Hi,Lu,} = {Hi,e,{K,)} = • {Hi,K,} - 0. 

So, the addition theorems (|2.6p could help us to classify algebraically superintegrable systems 
and vice versa. 

Remark 1 We can express hyperelliptic integrals i?y (|2.5ll2.6p via the elementary functions for 
the hyperelliptic curves Cj with genus =0 only [HH]. Namely, if 

degPj(A) >2, 

then variable loj (j2.5p are elliptic functions on p and q. 



Remark 2 According to [6J construction of the separated variables {q,p) for the L-subset of 
the Stackel systems is a pure computer problem. The software [6 may be easily improved in 
order to calculate additional integrals of motion ujk (|2.5|) too. 



Example 1 The (b) Drach system 

Let us consider the Stackel system defined by 



qt 92 
1 1 



Ui 2 = ±2a ■ 



Qt,2 



such that the separated relations (|2.3[) look like 



Pi 2 = Hiq^ ^ + H2±2a ^ . 

In variables /ij = qjpj and = q^ equations (|2.4[) became the canonical equations defining two 
Riemann surfaces: 



Ci_2 : = Pi,2(A) = -ffiA^ + H2X ±2aX-(3± 2-f. 

The action variables /1.2 = Hi 2 have the canonical brackets (jl.ip with the following angle 
variables (12.51): 



U!2 = 



A 



1 

VA(A) 



:dA 



92 A 



dA 



92 X 



:dA 



dA 



So, the Hamiltonian system with hamiltonian Hi has two independent integrals of motion H2 
(12:21) and 



1 



0J2 



where 



4^/i^^ 



P' - 



In piqi 



2VHiJ 



P292 



2\AHT 



dPi,2iX) 



dX 



2Hi qU+H2±2a. 



A=9? 



Using addition theorem 



e e = e 



or ln(a;i) + ln(a;2) = In (X1X2) 



(2.7) 



one gets 



W2 



1 



:ln 



2/Hi 



P2q2 



2VH1 
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So we can introduce algebraic integral of motion 

$i(/,c^) = e^^'^^ = ( p,q, + J%=] ( P2q2 ' 



4pip2gig2gi + P{P2 ^ piqiP-l + V2q2Pl ^ 



4Hi ' 

where K4 and are polynomials fourth and third order in momenta 

K3^2{pmP2+P2q2P[), K^^ PiP; + Apmp2q2Hi, (2. 

such that 

{m,K3} = {Hi,IU} = and {H2, K3} = K^. 
After the following change of variables 

(91-92)^ P1-P2 ((71+92)^ P1+P2 

X = -. > = : y = > Py = \ > 

4 qi - q2 4 qi + q2 

the Hamiltonian 

-Hi = PxPy + —==■ + 



xy {x-yY Jxy{x-yy 



coincides with the Hamiltonian for the one of the Drach systems (case (b) in [T^), whereas 
cubic polynomial K2 is the Drach integral of motion [T7] . 

Example 2 The (1) Drach system. 

Let us consider another Drach system (case (1) in [TT) with the Hamiltonian 



H, = 



PxPy + a {y - + I3x ^/^+7a; ^/^(y-px) 



Without lost of generality we can put p = —3. Substituting this Hamiltonian into the computer 
programm from [6] one gets the separated variables 

{qi-q2Y {qi+q2? 

and the corresponding separated relations 

P| = ^:) fe) =-4^91 T 8^27^3 + 4i7ig2^4V2/3gj +i?2, J = 1,2, (2.9) 



which give rise a pair of hypercUiptic curves at /ij — pj and Xj = qj . Using the Stackel matrix 

5* = 



Aqf Aql 



we can get the angle variable 



1 1 



1 dA 1 dA 



W2 = - 



27 /f^Aj 27 /fMA)' 

which is a sum of the incomplete elliptic integrals of the first kind. We do not know how to get 
algebraic integral of motion starting with this variable 0^2- 

However, we know that there is additional cubic integral of motion 

K3^2{P[p2 + P!2Pi). 
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which looks hke as the Drach integral (12. 8p . but in this case functions 

.2 d P,{\,) 



P' = (Ai + A2 



(Ai+A2)4 



■^1,2—91,; 



have not simple algebro-geometric explanation. Nevertheless, this polynomial integral i^s has 
to be a function on polynomials i?i,2 and elliptic hmction uj2. Since, the sum of elliptic integrals 

has to be some function on polynomials -ffi.2 and Kj,. 

It will be interesting to get this function explicitly and to understand why it exists. 

Example 3 The Henon-Heiles system 

Let us consider the Stackel system defined by 

"S* = J -1 ) ' ^^'^ " "^^'^ ' 

such that 

pi^2 ^Hi±H2 + aql^2 ■ 

The action variables 

Pi±p|_^M±M), = (2.10) 

2 2 2 2 ^ ' 

have the canonical brackets (jl.ip with the angle variables 1)2.5^ : 

1 / dA dA 



2 \] VaX' + H1+H2 J VaA^ + Hi - H2 

So, the Hamiltonian system with the hamiltonian Hi has two independent integrals of motion 
H2 (|2l0l) and 

F|^arcsin(zi), F^arcsin(2:i), k 

{Hi+H2Y'^ {H1-H2Y/ 



U2 = const I , ,1;, 777 M/6 I ' (2.11) 



where F{z, k) is the incomplete elliptic integral of the first kind and 



^1,2 



31/4 Y {Hi±H2y/^ 

After the following change of variables 



x±y 

91,2 = —^, Pl^2=Px±Py 



the Hamilton function Hi (|2.10|) 

2,2 a.^(a:^ + 3y^) 
Hi^Vx+Py ^ 

coincides with one of the Henon-Hciles hamiltonians. 

Of course, integral of motion uj2 (|2.1ip remains an elliptic function in any coordinates and 
we can not get single- valued integral of motion on the whole phase space. 
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3 The Sklyanin algebra 



In this section we study a class of finite-dimensional Liouville integrable systems described by 
the representations of the quadratic r-matrix Poisson algebra, or the Sklyanin algebra: 



{T{X),T{^l)}^[r{X-^i),TWT{^,)], 



(3.1) 



1 2 

Herer(A) = r(A) Id, r(^) = Id® r(/i) and r(A - ^) is a classical r-matrix [5]-[TU]. In the 
simplest case of the 4x4 rational r-matrix 



r(A-M) 



X — fi 



n. 



n 



/ 


1 











\ 








1 












1 










V 











1 





matrix T(A) depends polynomially on the parameter A 

" C(A) D(A) 

" aA" + AiA"-i + ... + A„ 
7A" + CiA"-i + ... + C„ 

The leading coefficients a,(3,j,d and 2n coefficients of the detT(A) 



(3.2) 



(3.3) 



/?A" 4-BiA 
(5 A" + DiX 



n-l 

ri-1 



d{X) = detr(A) = {aS - /37)A2" + QiA^"-! + • • • + Q 



2n • 



(3.4) 



are Casimirs of the bracket (|3.ip . Therefore, we have a 4n-dimensional space of the coefficients 
Ai,Bi, Ci and Di with 2n Casimir operators Qi, leaving us with n degrees of freedom. 

3.1 Open lattices 

For so-called open lattices independent Poisson involutive integrals of motion are given by the 
coefRcients of the entry A{X): 

A{X) = aA" + i/iA"-i + • • • Hn, {Hk,H^} - . (3.5) 

In the special action-angle representation ;19J, one has n pairs of the action-angle variables: 



A(/„ 



0, 



LUm = V ^lni?(/m), m=l,...,n. 



(3.6) 



having the standard Poisson brackets (|l.ip . 

Since the action variables Ik are zeroes of the polynomial 

n 

A{X) = aA" + i/iA"-i + • • • i/„ = a n (A - /„0, 

m— 1 

initial integrals of motion are elementary symmetric function on 

n 71 n 

Hi = — ^ H2 = ^ ^fc^m, • ■ • 5 = (—1)^ ^rr 

m— 1 k^rii m—1 



and the matrix J (|1.2p is equal to 

/ -1 
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In this case all the functions (|1.3p are functionally independent and one gets that open 
lattices related with the Sklyanin algebra are maximally superintegrable systems. 

We have to underline that functions ujk (|3.6p may be found without integration, i.e. using 
pure algebraic constructions. As for the (b) Drach system one gets 

{h,w,} = {hMB{Ij)} = 0, ^ {I,,B{I,)} = 0, 1^3. 

So, for any Hamiltonian li we have additional integrals of motion B{Ij), which are polynomials 
on momenta and the remaining action variables Ij. 



3.2 Periodic lattices 

For so-called periodic lattices integrals of motion are given by the coefficients of the trace of 
r(A): 

T{X)=tiT{X) = {a + 5)\" + HiX"-'^ + ---H„, {H„ H^} ^ . (3.7) 
According to [5] -[12], in this case the separated coordinates Ui,Vi defined by 

= 0, Vi = A{ui) , i = l,...,n, 

satisfy to the following separated relations 

v^ + det T {u^)v~^ = A{u^) + D{u^) ^t{u^), i^l,...,n. (3.8) 

The corresponding symplectic form in (u, u)-variables is written as 

n 

n^ri"^d log(wj) A duj . (3.9) 

Let us take the action coordinates Im = H,n. (|3.7p instead of w- variables [TOl [11] . From (|3.8p 
one gets 

Vj = A{uj) ^ i {t{u.^ + y^T(uj) -4d(Mj)^ 
and one easily find expression for the symplectic form in (/, u)-variables 

j=i fc=i v^m) 

where P(A) ~ r^(A) — 4d{X). The equations of motion 

A - Uk 



{TiX),u,}^rjJpiu,) n 



are linearized by the Abel transformation [2l |4] and the corresponding angle variables look like 
u;j=vY, -7^'^^' degP(A)=2n. (3.10) 



k=l 



In this case matrix J ()1.2p is identity matrix and integrals of motion 0^ = ujk (|3.10p are elliptic 
functions and we could not get single-valued additional integral of motion on the whole phase 
space. 

Recall, that quadratic r-matrix algebra (|3.ip is related with a set of integrable systems that 
includes XXX Heisenberg magnet, the Toda lattices, the discrete self-trapping model and the 
Goryachev-Chaplygin gyrostat [19] . 
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Example 4 The Heisenberg magnet 

Our main model turns out to be an n-site Heisenberg magnet, which is an integrable lattice of 
n sl(2) spins with nearest neighbor interaction. In the lattice representation the matrix T(A) 
(|4.2|) acquires the following form: 



T(A) = Li(A - ci) L2{X ~ C2) • • • L„(A - c„) , (3.11) 

with 

Here the local variables sf , i — 1, . . . ,n, are generators of n copies of the sl(2) Poisson algebra: 

{4, sf} = ±sf, {s+ sj} = 24 . (3.12) 
and Cm are arbitrary numbers. If n = 2 one gets 

A{X) = X^+ HiX + H2 = {X- h){X - I2), 

where 

Hi = -{4 + sl+Cl+ C2), H2 = slsl + S+S2 + Cisl + C2sl + ClC2 , 

and 

1 



ha = -^ [Hi±^H(-4H2 
The second polynomial 

B{X) = (sf + s^)X - sjs^ + s^sl - cisf - C2sf 
defines the angle variables 

wi,2 = lnB(/i_2) . 

The functions (j)i^2 p.3p 

hoJl - I2OJ2 . UJl ~ UJ2 



'■>! — 



h-h ' h~h 



obey to equations {Hi^cfyj} = Sij. So, the Hamiltonian system with the hamiltonian H2 has 
two independent integrals of motion Hi and 

$ = -(/i - I2) 01 = /iwi - /2W2 = h lnB(/i) - /2 lnB(/2) . 

The similar additional integral for the periodic lattice is some combination of elliptic functions 
even in the simplest case n — 2. 

Example 5 The Toda lattice 

The Toda lattices appear as a specialization of our basic model when the parameters are fixed 
as follows: 

/? = 7 = (5 = and det r(A) = 1. (3.13) 

We also put a — 1 and rj — —1. In the lattice representation, the monodromy matrix T (j4.2p 
acquires the form 

r(A) = ii(A) • • • L„_i(A) i„(A) , L,^(^ ^ . (3.14) 

Here pi, qi are Darboux variables {qi,pj} = Sij. If n = 2 one gets 

.4(A) = A2-(pi+p2)A+piP2-e«^-«^ = (A-/i)(A-/2), 
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and 

Let us consider natural Hamiltonian systems with integrals of motion 

Hi = -Hi=pi+p2, ff^ ^ Pi +P2 ^ cgi-g2 



2 



2- 



For these integrals we have 



1 1 \ ~ I1UJ2 - h^i ~ UJ1-UJ2 

J = I T T I , and cpi = 



h I2 J ' Ii ~ I2 ' " Ii — h 

So, the Hamiltonian system with the hamiltonian H2 has two independent integrals of motion 
Hi and 

$1 = -(/i -l2)4>i =hoJ2-hoJi = -{hlnBih) - hl-aBih)) . 

Namely this " generalized angular momentum" has been found in [1] [3j . In contrast with the 
Drach system we can not directly apply additional theorem in this case. 
At the periodic case one gets 

h = ~{pi + P2), h = P1P2 - e"'-"' - e«=-«S 

and the Hamiltonian system with the hamiltonian H = I2 has two independent integrals Ii and 

gi + 92 , hw2 , Z"^' dA 

(pi = 1 — , where uj2 - 



V(A2+/iA + /2)2-4 

As for Henon-Heiles system this additional integral of motion may be rewritten as a combination 
of elliptic functions on pi and Ii 2- 

The global action-angle variables for the periodic Toda lattices are discussed in ^5]. Nev- 
ertheless we suppose that periodic Toda lattices are locally maximally superintegrable systems 
only. 



4 Reflection equation algebra 

In this section we study a class of finite-dimensional Liouville integrable systems described by 
the representations of the reflection equation algebra: 

{r(A),r(A.)} = [r(A-M),T(A)r(M)] (4.1) 

+ r(A)r(A + ^)f(/i) -T{fi)r{\ + n)T{\) . 

In the simplest case of the 4x4 rational r-matrix (j3.2p matrix T(A) depends polynomially on 
the parameter A 

f A{X) B{\) \ , ( 2n+l 2n + l\ 

Coefficients of the entries 

^(A) = a A2"+i + A2n A'" + ^2„-i A2"-i . . . + ^0, 

B{\) = A2«+i + 6„A2"-i + 6„_iA2"-3 _ ^ g^^^ (4 3) 

C(A) = C„A2"-i + . . . + CaA-^ + CiA, 



9 



are generators of the quadratic Poisson algebra (|3.1|) . The leading coefficient a and 2?! + 1 
coefficients of the detT(A) 

d{X) = det T(A) - Q2„A^" + Qa^-iA^''-^ + ■ ■ ■ + Qo ■ (4.4) 

are Casimirs of the bracket (|3.1[) . Therefore, we have a 4n + 1-dimensional space of the coeffi- 
cients 

■^Ot ■ ■ ■ 1 ^2n: ^1: ■ ■ ■ , ^n, Ci,..., Cn (4-5) 

with 2n + 1 Casimir operators Qi, leaving us with n degrees of freedom. 
4.1 Open lattices 

For so-called open lattices integrals of motion are given by the coefficients of the entry B{X): 

S(A) = A2"+i+i/iA2"-i+if2A'"-3 ..-fff„A, {Hk,H,-a} = 0. (4.6) 
In the special action-angle representation ^20], one has n pairs of the action-angle variables: 



B{±lm) = 0, uJm = V ^lnA{Im), m = 1, . . . , n, 

having the standard Poisson brackets (jl.ip . 

As above the action variables Ik are zeroes of the polynomial 



(4.7) 



B(A) = A2"+i + i/iA^"-! + i?2A2"-3 . . . + H„X = A [| (A^ - /, 

m— 1 

integrals of motion Hm are elementary symmetric function on 

71 n n 



m— 1 k^m 



2 

m 1 



and the matrix J (|1.2p is equal to 

/ -^1 



J = 2 



m^l 

(-i)"^i n 



-In 

I ^ l'^ 

m=^n 



i-^Th n Irn 



As above all the functions </ij (jl.3p are functionally independent and one gets that open lattices 
related with the reflection equation algebra are maximally superintegrable systems. 

We have to underline that functions LOk (|4.7p and, therefore, integrals of motion (pk may 
be found without integration, i.e. using pure algebraic constructions. As above any action 
variables Li have additional polynomial integrals of motion A{Ij) 

{h,w,} = {hMA{Ij)} = 0, ^ {I,,AiIj)}^0, i^j, 

but for the integrals Hk one gets some combinations of logarithms and Ij only. 
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4.2 Periodic lattices 



The theory of periodic lattices is based on the following construction of commutative subalgebras 
[9l I12j . Let us introduce the boundary matrix 

IC(X) = ( "f^J J,, ) (4.8) 



c(A) d(A) 

whose entries a(A), c?(A) are polynomials with numerical coefficients and entry c(A) is arbitrary 
polynomial on A. If the polynomial 

t(A) = tr /C(A)T(A) = a{X)A{X) + c{X)B{X) + d{X)A{-X) 

has n independent dynamical coefhcients Hi , . . . , Hn only, then 

{r(A),r(Ai)} =0, ^ {H,,H,} = 0, t,j^l,...,n. 

These Poisson involutive integrals of motion Hi define the Liouville integrable systems or the 
periodic lattices related with the reflection equation algebra. 

In the periodic case the n pairs of the former action-angle variables (|4.7p for the open lattices 

S(±u„i) = 0, Vm^A{um), m=l,...,n, (4.9) 

are the simple separated variables, which satisfy to the following separated relations 

a{uj)vj + d{uj) detT(uj) vj^ = T{uj) . (4-10) 

The corresponding symplectic form in (u, w)-variables has the form p.9p . Therefore, we can ap- 
ply the above consideration of periodic lattices related with the Sklyanin algebra to construction 
of the angle variables (|3.10p in this case too. 

As above the action variables Im coincide with integrals of motion Hm and the matrix J 
(|1.2p is identity matrix. Additional integrals of motion (pk coincide with the angle variables LUm 
(|3.10p . which are integrals of Abelian differentials on the hyperelliptic curve = P{X) defined 
by (I4.10p . i.e. they are elliptic functions. So, all the periodic lattices related with the reflection 
equation algebra are formally maximally superintegrable systems, i.e. locally superintegrable 
only. 

Recall, that reflection equation algebra (|4.ip is related with a set of integrable systems that 
includes the generalized Toda lattices, the Kowalevski top and XXX Heisenberg magnet with 
boundary conditions pU] , 

Example 6 Open Toda lattice 

According to [51 [H] the 2x2 Lax matrix for the generalized open Toda lattice 



T{X) =[11 W)] IC-W ( n I ' (4.11) 

where Li is given by (|3.14p and 



^fe=i / \fc=i 



1^ (^\_( 2a2A2 - ioiA -I- oo (4a2e«" - 1)A 
' ^ \ 2a2A2 + zaiA + ao 

satisfies to the reflection equation algebra at a = and r/ = 1. Here pt^qi are dynamical 
variables and are parameters. 

All the open Toda lattices associated with classical root systems "^n, and Sin are 

isomorphic to each other [T^ and, therefore, we consider ^„ root system only. In this case 
a2 = ai — and at n = 2 one gets 

B^X^ + HiX^ + H2X = A(A - /2)(A - li) = 

= A^ - {pI +pI- 2e9i-«^ - 2aoe«^)A3 + {{pip2 - e'^'-^'f - 2ao{plei' + e«0)A, 
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such that 



^1— ^^1^^2' ^2— -'1^2: (pl — - rt/T2 72\~' 



So, the Hamiltonian system with the hamiltonian Hi has two independent integrals of motion 
H2 and 

$2 = 2hh{l! - Il)4>2 = huJ2 - /2Cc'l = -(/l hl^(/2) - /2 ln^(/2)), 

where 

A = -e-«i (a* - piA^ - {pI + 2aoe«^ + + (pip^ _ pae^^-"?^ + 2aoPie«^ )a) -oq . 

Such as matrix J (|1.2|) has the special form, the additional integral of motion $2 has the form 
of the angular momentum in (/, w)-variables (see ^I7j and |3j). However at n > 2 we will have 
completely another picture. 



5 Conclusion 

We discuss some constructions of the single valued integrals of motion on the whole phase space 
by using multi-valued action-angle variables. 

For the Stackel systems we use addition theorem for the construction of the polynomial 
additional integrals of motion starting with zero-genus hyperelliptic curves. The construction 
of the polynomial additional integrals of motion is an algebraic procedure for the open lattices 
related with the Sklyanin algebra or with the reflection equation algebra. On the other hand 
for the periodic lattices one gets the sums of integrals of first kind Abelian differentials on the 
hyperelliptic curves, which is locally defined only. 

It will be interesting to consider inverse problem and try to find algorithm of construction 
action-angle variables starting with known additional polynomial integrals of motion. 

The research was partially supported by the RFBR grant 06-01-00140. 
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